Spin Polarization of Gapped Dirac Surface States Near the Topological Phase 

Transition in TlB^Si^Se^ 
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We performed systematic spin- and angle- resolved photoemission spectroscopy of TlBi(Si_ I Se :c )2 
which undergoes a topological phase transition at x ~ 0.5. In TlBiSe2 [x = 1.0), we revealed a 
helical spin texture of Dirac-cone surface states with an intrinsic in-plane spin polarization of ~ 0.8. 
The spin polarization still survives in the gapped surface states at x > 0.5, although it gradually 
weakens upon approaching x = 0.5 and vanishes in the non-topological phase. No evidence for 
the out-of-plane spin polarization was found irrespective of x and momentum. The present results 
unambiguously indicate the topological origin of the gapped Dirac surface states, and also impose 
a constraint on models to explain the origin of mass acquisition of Dirac fermions. 

PACS numbers: 73.20.-r, 71.20.-b, 75.70.Tj, 79.60.-i 



Three-dimensional topological insulators (TIs) exhibit 
a novel quantum state with metallic topological surface 
state (SS) which disperses across the bulk band gap gen- 
erated by a strong spin-orbit coupling. Specifically, the 
topological SS is characterized by a linearly dispersing 
Dirac-cone energy band with the helical spin texture [J- 
4], which hosts massless Dirac fermions protected by 
time-reveal symmetry (TRS). This peculiar SS of TIs 
provides a platform for fascinating quantum phenomena 
such as the robustness against nonmagnetic impurities / 
disorder [B], Q and the emergence of Majorana fermions 
0. Breaking the TRS by a magnetic field or a mag- 
netic order lifts the band degeneracy at the Dirac point 
and causes an energy gap called a Dirac gap in the topo- 
logical SS, turning the massless Dirac fermions into a 
massive state. Theoretically, realization of the massive 
state is a prerequisite to novel topological phenomena 
such as the topological magnetoelectric effect and half- 
integer quantum Hall effect [1, Q, but experimentally, 
the role of TRS-breaking magnetic impurities in trigger- 
ing the Dirac-g ap ope ning in the topological SS is under 



intense debate lO- 1 



Recently, it was found by angle-resolved photoemission 
spectroscopy (ARPES) measurements that the thallium- 
based ternary chalcogenide solid-solution TIBi(Si_ x Se x )2 
goes through a quantum phase transition (QPT) from 
the topolo gica l [x > 0.5) to the non-topological (x < 



0.5) phase [lj,ll5[. Intriguingly, the Dirac-cone topolog- 
ical SS in TlBiSc 2 (x = 1.0) starts to show a Dirac gap 
upon partially replacing Se with S, and the gap gradually 
grows upon increasing the S content until the QPT at x 
= 0.5 eliminates the SS jl5| . Since there is no evidence 
for magnetic impurities in the measured crystals, it has 
been puzzling how the Dirac fermions in TlB^Si-xSe^ 
acquire the mass without explicitly breaking the TRS. 
To gain insights into the microscopic mechanism of this 



intriguing mass acquisition, it is of particular impor- 
tance to establish the spin structure of the gapped SS in 
TlBi(Si_ a; Se 2 ;)2, since the topological or non-topological 
nature of the SS as well as the origin of the Dirac gap is 
intimately related to the characteristics of the spin tex- 
ture. 

In this Letter, we report our spin-resolved ARPES re- 
sults of TlBi(Si_ x Se x )2 with various x values [x = 1.0, 
0.9, 0.8, 0.6, and 0.0), and present the evolution of the 
spin texture as a function of x. We demonstrate that the 
helical spin texture persists even in the gapped SS ob- 
served at x > 0.5, while the spin polarization vanishes in 
the ARPES spectrum in the non-topological phase at x 
< 0.5 in accordance with the absence of the SS. In addi- 
tion, we found unusual reduction of the spin polarization 
in the SS near the QPT. We discuss the implications of 
the present experimental findings on the mechanism of 
the Dirac gap opening. 

High-quality single crystals of TlBi(Si_ a; Se a ;)2 were 
grown by a modified Bridgman method. Details of 
the sample preparations were described elsewhere [lij ]. 
ARPES measurements were performed with a spin- 
resolved photoemission spectrometer based on MBS-Al 
analyzer at Tohoku University [l6(. We used one of the 
Xe I lines (hv = 8.437 eV) to excite photoelectrons. Sam- 
ples were cleaved in situ along the (111) crystal plane in 
an ultrahigh vacuum of 5xl0 -11 Torr. The energy reso- 
lution for the spin-resolved and regular ARPES measure- 
ments was set at 40 and 6 meV, respectively. The sample 
temperature was kept at 30 K during the measurements. 
We used the Sherman function value of 0.07 to obtain 
spin- resolved ARPES data fl6|. 

First we demonstrate the spin structure of TlBiSe2 (x 
= 1.0) in which the topological origin of the Dirac-cone 



SS is well established by previous studies [17H19j . Figure 



1(a) shows the ARPES intensity plot at E F of TlBiSc 2 
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FIG. 1: (Color online) (a) ARPES intensity at Ep around the 
f point for TlBiSe2 plotted as a function of two-dimensional 
wave vector k measured with the Xe I resonance line (hv = 
8.437 eV) at T = 30 K. The ARPES intensity is integrated 
within ±20 meV with respect to Ef. The definition of an 
in-plane up-spin vector is also indicated by an arrow, (b) 
Near-£ F ARPES intensity along the YK cut for TlBiSe 2 . (c, 
d) Spin-resolved EDCs for the k points (A-E and F-J) whose 
positions in the 2D Brillouin zone are indicated by open circles 
in (a). 



around the f point as a function of in-plane wave vector 
k x and k y . We clearly observe circular Fermi surface (FS) 
with hexagonal deformation, which originates from the 
upper branch of the Dirac-cone SS as displayed in Fig. 
1(b). We also find a finite intensity distribution around 
the r point as indicated by white arrow in Fig. 1(b). 
This signal is attributed to the bulk conduction band 
(CB) which is slightly occupied due to electron doping 
caused by Se vacancies in the crystals 
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Figures 1(c) and (d) show spin- resolved energy distri- 
bution curves (EDCs) of TlBiSe 2 for the in-plane spin 
component measured at various k points along the TM 
cut as indicated in Fig. 1(a). We define the spin po- 
larization vector perpendicular to the measured k [see, 
e.g., thick red arrow in Fig. 1(a)], where the "up spin" 
points to the clockwise direction. One can immediately 
recognize that the up-spin EDC is always more promi- 
nent than the down-spin counterpart in all the measured 
k region, consistent with the spin helical texture of the 
Dirac-cone SS 0, E3] • Upon closer look, one finds that 
the shape of the up-spin EDC depends on the k point. 
For instance, near the T point (point E and J) the peak 
in the up-spin EDC (marked by circle) is located at ~ 
0.3 eV, but it approaches Ey on moving away from the 
f point (from E to A, or from J to F), reflecting an elec- 
tronlikc dispersion of the upper branch of the Dirac cone 
(note that the dispersive feature is less clear than the reg- 
ular ARPES data because of the wider k window and the 
lower energy resolution of the spin-resolved mode). On 




(c) 



0.4 0.2 E f 0.4 0.2 E f 04 0.2 E f 0.4 0.2 E, 0.4 

Binding Energy (eV r Binding Energy (eV) Binding Energy {eV] Binding Energy (eV) Binding 



0.2 E f 
Energy (eV) 



S 0.0 
o 

0. 

■0.3 




' o!4 ' 0.2 ' 4 04 ' 6.2' 5 0i4 0.2 E, 0i4 0.2 E, 0.4 0.2 ^ 

Binding Energy (eVJ Binding Energy (eV) Binding Energy reV) Binding Energy feV! Binding Energy (eV) 

FIG. 2: (Color online) (a) ^-dependence of near-Ep ARPES 
intensity around the F point for TlBi(Si_ x Se a: )2. Dashed lines 
are guides to the eyes to trace the band dispersion. The ver- 
tical red (thick) line represents the k point where the spin- 
resolved EDCs were obtained, (b, c) Corresponding spin- 
resolved EDCs for the specified k point and their spin polar- 
ization Pip, respectively. In (c), the averaged Pip value in the 
energy window of ±0.1 eV with respect to the energy of the 
SS (Ess) is indicated by a number and also by a dashed line. 
The experimental uncertainty (error bar) in the averaged Pip 
value is ±0.05. 



the other hand, the peak in the down-spin EDC is always 
pinned at Ep. This is because it is governed by the bulk 
CB which contributes equally to the up- and down-spin 
EDCs around the f point (as evidenced by the small dif- 
ference in the near-Ep spectral weight between the up- 
and down-spin data at points D, E, I, and J) as a result 
of the spin-degenerated character. 

Figure 2(a) shows the near-_Ep ARPES intensity plot 
of TlBi(Si_ a Se 2; )2 measured along the TM direction (k x ) 
for various x values (x = 1.0-0.6 in the topological phase 
and x = 0.0 in the non-topological phase). The "X"- 
shaped topological SS transforms into the gapped SS in 
x = 0.9-0.6, as one can see from the suppression of the 
ARPES intensity at the Dirac point and the parabolic- 
like dispersion of the upper branch of the Dirac cone. The 
characteristics of the gapped SS and its intrinsic nature 
are demonstrated in Supplemental Materials, where it 
is explicitly shown that the Dirac gap opening in our 
data is not due to a spectral-weight suppression at the 
Dirac point, momentum broadening, bad cleaves, nor a 
random disorder effect. At x = 0.0 (TlBiS 2 ), the SS 
completely disappears and only the bulk CB is observed 
at Ep, reflecting the non-topological nature [HI]. 

To elucidate the evolution of the spin structure of the 
SS, we have performed spin-resolved ARPES measure- 
ments by focusing on a particular k point. The k point 
was chosen to be located midway between the Fermi vec- 
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FIG. 3: (Color online) (a) Representative results of the two- 
component fitting to the EDC obtained by the regular ARPES 
experiment to estimate the bulk and surface spectral weights 
as indicated by blue and red curves, respectively. The plot- 
ted EDCs are generated by summing the raw EDCs within 
the k window of the spin-resolved measurement, to take into 
account the k broadening effect. In the fitting, we assumed a 
Lorentzian peak for the surface, bulk VB, and bulk CB, mul- 
tiplied by the Fermi-Dirac distribution function. A constant 
term was subtracted from the peaks for bulk VB and CB to 
reproduce the bulk band gap. (b) Ratio of the surface spectral 
weight to the total spectral weight, Nss/N to tai (where -/V tota i 
= Nss + A?buik), at Ess plotted against x as estimated by the 
fitting of EDCs [see Fig. 3(a)]. (c) P IP divided by iVss/iVtotai, 
which represents the intrinsic spin polarization Pj nt of the SS 
(circles). Pip is also plotted with triangles. 



tor (fcp) and the T point along the f M direction [see 
vertical red (thick) line in Fig. 2(a)] so that the binding 
energy (£"b) of the surface band (Ess ~ 0.25 eV) is apart 
from that of the bulk CB. As can be seen in Fig. 2(b), 
the shape of the spin-resolved EDCs for the in-plane com- 
ponent is similar between x — 1.0 and 0.9. In addition, 
the spin polarization (Pip) [Fig. 2(c)] exhibits a simi- 
lar value of ~ 0.3, suggesting that the gapped SS in x 
= 0.9 has a helical spin texture as in x — 1.0. Another 
notable aspect is that the in-plane spin polarization is 
gradually reduced upon approaching the QPT at x ~ 
0.5, as evidenced by the smaller difference between the 
up- and down-spin EDCs for the samples with lower x 
values. The Pip value in x — 0.6 is ~ 0.15 (half of that 
at x = 1.0), and it vanishes in the non-topological phase 
(x = 0.0). 

To discuss more quantitatively the absolute spin po- 
larization of the SS, it is necessary to remove the finite 
contribution from the spin-degenerated bulk band. For 
this purpose, we have fit the EDCs in the regular ARPES 
measurement by incorporating both the surface and bulk 
components [see Fig. 3(a)], and estimated the spectral 
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FIG. 4: (Color online) (a) Spin-resolved EDCs and corre- 
sponding spin polarization for the out-of-plane component in 
the topological phase (x — 1.0, 0.8, and 0.6) measured at the 
k point along the TM direction indicated by the red (thick) 
line in Fig. 2(a). The EDCs along the FK direction are also 
plotted for x_= 1.0. (b) Near-£ F EDCs for x = 0.9. Spec- 
trum at the F point is highlighted by thick black curve, (c) 
Spin-resolved EDCs and corresponding spin polarization for 
the out-of-plane component measured at the F point for x — 
0.9. (d) Schematic picture of the evolution of the topological 
SS (red) and bulk bands (blue) as a function of x. The spin 
polarization vector is indicated by black arrows. 



weight of the SS (A^ss) with respect to that of the total 
weight (Atotai = N S s + -Nbulk)- The intrinsic in-plane 
spin polarization of the SS (Pint) corresponds to the Pip 
value divided by Nss/N to tai- (For the detailed procedure 
to estimate the spin polarization and corresponding error 
bars, see Supplemental Materials). Figures 3(b) and (c) 
display the x- dependence of Nss/N to tai and P mt (Pip is 
also plotted for reference), respectively. As can be seen 
in Fig. 3(c), the estimated P n t value, which reaches ~ 
0.8 for x = 1.0, is gradually reduced with increasing x to 
become ~ 0.5 at x = 0.6, and finally vanishes in the non- 
topological phase at x = 0.0. The origin of the reduction 
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in the in-plane spin polarization toward the QPT in the 
gapped SS is not clear, but one can infer that it is not 
due to contamination from bulk conduction bands or the 
degenerate surface two-dimensional electron gas (2DEG) 
[20l | . since the effect of the former is already discounted 
in our analysis and the latter is obviously not formed in 
our surface situation. The less-than-unity value of P; n t at 
x = 1.0 is likely to be a consequence of the spin-orbit en- 
tanglement as proposed from first-principle calculations 

The next important issue to clarify is the out-of-plane 
(z) spin component. Figure 4(a) displays the spin- 
resolved EDCs and corresponding spin polarizations for 
the z-spin component measured at the same k point as 
in Fig. 2 (along the TM direction) for the samples in the 
topological phase (x = 1.0, 0.8, and 0.6). Spin-resolved 
EDCs for the k point along the f K direction are also 
displayed for x = 1.0; obviously, there is no apparent z- 
spin component in the TM nor the TK directions. It was 
previously found @, [22|, [23| that in Bi 2 Te3 a finite z-spin 
texture emerges along the TK direction due to strong 
hexagonal warping of the SS. However, the hexagonal 
warping effect is much weaker in TlBi(Si_ 2; Se 2: )2 than 
in Bi2Te3, resulting in no or negligible z-spin component 
The present data reconfirms this previous result for 
the whole topological phase, where the estimated z-axis 
spin polarization remains negligible (-0.02 ~ +0.01). 

While the above result demonstrates no z-spin polar- 
ization near the &f point, a natural question arises as 
to whether any z-spin component exists exactly at the 
r point. To address this issue, first we carefully exam- 
ined the EDCs for x = 0.9 [Fig. 4(b)] and determined the 
spectrum which exactly corresponds to the T point (thick 
line) , and then performed the spin- resolved ARPES mea- 
surement at this k point with a higher energy resolution 
(15 meV). As shown in Fig. 4(c), we have succeeded in 
resolving two well-separated peaks in the spin-resolved 
EDCs attributed to the upper and lower branches of the 
gapped SS, respectively. Notably, there is essentially no 
difference between the up- and down-spin EDCs, indicat- 
ing null z-spin polarization. These results suggest that 
the z-spin component is absent irrespective of x and k. 
Based on these observations, we depict in Fig. 4(d) the 
schematic picture of the spin-texture evolution for three 
different phases (massless Dirac, massive Dirac, and non- 
topological phases) in TIBi^i^Se^- 

Now we discuss the implications of the observed char- 
acteristics of the spin polarization in relation to the 
Dirac-gap opening. First, the present result directly ad- 
dresses the question whether the gapped SS has a topo- 
logical origin. One may argue that the topological phase 
transition should have occurred between the massless (x 
— 1.0) and the massive (x = 0.9) phase, and that the ob- 
served gapped SS corresponds simply to a pair of trivial, 
spin-degenerate SSs which happened to present a gap- like 
feature at the f point. In this respect, such a possibility 



is immediately ruled out by the observation of the spin 
non-degenerate nature of the gapped SS. 

Another possible origin of the gapped SS is the hy- 
bridization of two topological SSs. In ultrathin films of 
Bi2Se3, it was reported that a Dirac gap opens when the 
electron wavefunctions of the top and bottom surfaces 
start to hybridize 3, 2f|. In this case, the hybridized 
state is spin degenerate and loses the original helical po- 
larization [2^]. While the thick nature of our samples 
(~ 100 jum) already spoke against this possibility |15| . 
the spin non-degeneracy observed here conclusively rules 
out the hybridization scenario. Moreover, a recent the- 
ory for ultrathin films in the hybridized state predicted 
a finite z-spin component [25| . which was not observed 
in the present experiment. It is worth noting that re- 
cent ARPES studies on Bi2Se3 ultrathin films 27J and 
bulk heterostructures involving Bi2Se3 units [28[ have re- 
ported a second Dirac cone which originates from the 
backside surface (interface). If the hybridization of the 
wave function between the Dirac cone at the topmost 
surface and the second Dirac cone at a backside surface 
(which may accidentally form due to bad cleaves) were 
the source of the observed phenomena, the separation be- 
tween the topmost surface and the accidentally-formed 
backside surface must become systematically closer with 
decreasing x, which is extremely unlikely. 

An important question is whether the TRS is broken in 
the massive Dirac phase. When the topological SS feels 
a magnetic field or magnetization directed perpendicu- 
lar to the surface (along the z-axis), it obtains a Dirac 
gap and the spin vector at the T point is aligned along 
the z-axis in opposite directions in the lower and upper 
branches of the gapped SS Q. In this respect, our finding 
that there is no z-spin component in the gapped SS does 
not give any support to this scenario. Nevertheless, one 
should keep in mind that if there are magnetic domains 
with antiparallel polarizations, the overall z-spin compo- 
nent can be canceled. Hence, it is difficult to draw any 
firm conclusion about the TRS breaking from the present 
experiment. 

A viable possibility for the origin of the Dirac gap is 
the bulk-surface coupling. A recent theory suggested that 
the Dirac-cone feature can be destroyed by strong poten- 
tial impurities which lead to bulk-assisted scattering to 
allow virtual spin- flip excitations (2i|. The proximity of 
the bulk band to the Dirac point near the QPT and 
the disorder due to random replacement of S with Se 
may be in line with this scenario. Another possibility 
could be critical fluctuations associated with the QPT as 
mentioned in Ref . [15| , and the reduction of the net spin 
polarization toward the QPT may be in line with this sce- 
nario; however, no evidence for such critical fluctuations 
has been obtained so far. While a comprehensive answer 
to the origin of the Dirac gap has not yet been obtained, 
the present observation that established the nature of 
the spin texture and its evolution in the gapped phase 
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provides an important step toward microscopic under- 
standing of the topological QPT and the unconventional 
mass acquisition of Dirac fermions in TIs. 

In summary, we have reported our spin-resolved 
ARPES experiments on TlBi(Si_a;Se s )2 to determine the 
spin polarization as a function of x. Besides firming up 
the intrinsic nature of the Dirac gap observed for 0.6 
< x < 0.9, we found that the helical spin texture sur- 
vives in the massive Dirac phase, which clearly points to 
the topological nature of the gapped SS. While the exact 
mechanism of the Dirac gap opening is yet to be eluci- 
dated, the direct connection between the emergence of 
the helical spin texture and the topological QPT, as well 
as the reduction of the net spin polarization toward to 
the QPT, gives an important clue. 
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help in sample preparations. This work was supported by 
JSPS (NEXT Program and KAKENHI 23224010), JST- 
CREST, MEXT of Japan (Innovative Area "Topological 
Quantum Phenomena"), AFOSR (AOARD 104103 and 
124038), and KEK-PF (Proposal number: 2012S2-001). 
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SI. Characteristics of the Dirac gap 



To elucidate the characteristics of the gapped surface state (SS) in more detail, we show in Fig. SI the energy 
distribution curves (EDCs) at the T point for TlBi(Si_ a; Se a ;)2 at each selenium composition x in the topological phase. 
As one can see in this plot, the EDC at the f point for x — 1.0 is composed of a single peak at the Dirac point (DP) 
which is satisfactorily fitted by a single Lorentzian peak. On the other hand, we observe two well-separated peaks for 
x = 0.9-0.6 (though the peaks for x = 0.6 are broad), which are very well reproduced by two Lorentzian peaks. If the 
energy gap of the SS (the Dirac gap) were caused by a spectral- weight suppression at the DP as proposed in Ref.[30|, 
there is no reason that the experimental data presenting a two-peak feature (particularly the data points between the 
two peaks) are fitted well by a pair of Lorentzians. Our result thus suggests that the observed two-peak structure 
reflects an intrinsic spectral function of the SS with an energy gap at the DP. 

We note that the data for the gapped phase displayed in Fig. SI were measured independently from our previous 
published work [l5| using different samples, and yet, the Dirac gap shows a quantitatively similar size indicative of 
its reproducibility. Specifically, the gap size in our previous experiment 15], as estimated from the peak position 
of the second-derivative of the EDCs, were 50 ± 10, 70 ± 10, 130 ± 20 meV for x — 0.9, 0.8, and 0.6, respectively, 
and those in the present study as estimated from the energy separation of the peak position of two Lorentzians are 
57 ± 2, 73 ± 2, 142 ± 5 meV, respectively, which agree well within the experimental uncertainties. In fact, as shown in 
Fig. S2 which plots the cleaving dependence of the band dispersion for x = 0.9, the obtained gap feature is robust: 
the gap always exhibits the same value (~ 55 meV) despite the finite variations in the overall intensity distribution for 
different cleaves. This result is obviously incompatible with the recent proposal [3fJ that the "gap" can be caused by 
bad cleaves and rough surfaces, which would lead to sizable variation in the gap value for different cleaving conditions. 
It should also be noted that if the gap is caused by a momentum broadening, the spectral linewidth for the gapped 
sample should be broader than that for the samples showing no gap. However, as demonstrated in Fig. S3, the peak 
width of the momentum distribution curves (MDCs) for x = 0.9 is comparable to or even narrower than that for x — 
1.0, indicating that the momentum-broadening effect is unlikely to be the source of the gap opening. 

One may argue that a "gap" could be caused by strong disorder and resultant smearing of the spectral weight 
at the DP. If this is the case, it is naturally expected that a similar gap appears in other solid-solution systems of 
topological insulators. In this respect, it is useful to compare the present result with Bi2- a; Sb 2 ;Te3_j ) Sey (BSTS) [3~0 | 
and Pb(Bii_ 2 ,Sb 2 ;)Te4 (Pbl24) [32j where similar or even stronger random disorder is introduced. As one can see 
in Fig. S4, we found no evidence for a Dirac gap in both BSTS and Pbl24 systems regardless of the ordered or 
disordered compositions, since the upper and lower branches of the Dirac cone are obviously degenerate exactly at 
the r point. Hence, in our experiments the Dirac gap has so far been only observed in the TlBi(Si_ 2; Se a ;)2 system in 
proximity to a topological QPT. 



S2. Procedure to estimate the spin polarization and its accuracy 



In the analysis of the spin-resolved ARPES data, we first obtained the scattering asymmetry parameter A = 
(Al — Ar)/(Al + Ar) at each binding energy where Al and Ar are the number of electrons detected by the 
left-hand-side and right-hand-side channeltrons in the Mott spin polarimeter, respectively. Next we obtained the 
in-plane spin polarization Pip through the formula Pip = A/S c s where S e g is the effective Sherman function (S e ff 
= 0.07; Ref.0). The statistical error in Pj P (APjp) is directly linked to that in A(AA) through APrp = AA/S eS . 
Since AA ss AA/A(A = Al + Ar), we accumulated all the spectra in Figs. 2(b) and (c) until the same error-bar 
range (±5%) is achieved for the Prp value when it is averaged over the energy region of ±0.1 eV centered at the 
surface-band energy. We plot this averaged P\p value as a dashed line in Fig. 2(c) (with its actual value indicated), 
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and also by triangles (with error bars of ±5%) in Fig. 3(c). The error bars in the intrinsic spin polarization Pj nt 
[-Pint = -Pip / (-^ss/^totai)] were evaluated from that in Pip and Ngs / N to ta,i by taking into account the error propagation. 
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FIG. SI: (a) EDCs at the f point for various x values measured in TlB^Si-nSe^. Black bar indicates the peak position of 
the Dirac-cone surface state estimated by the numerical fitting, (b) Results of fittings to the EDCs shown in (a). We assumed 
the Lorentzian peak(s) for the surface state (blue and yellow solid curves) and a polynomial function of up to the 4th order for 
the background (dashed curve). The size of the Dirac gap, as estimated from the energy difference between the peak positions 
of two Lorentzians, is also indicated for x = 0.9, 0.8 and 0.6. 
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FIG. S2: (a) ARPES intensity plots around the DP of TlBi(So.iSen.9)2 (x = 0.9) for five different cleaved surfaces, and (b) the 
corresponding EDCs. Dashed curves in (b) are guides to the eyes to trace the surface-band dispersion. The size of the Dirac 
gap, as estimated with the same procedure as in Fig. SI, is indicated by black solid lines. 
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FIG. S3: Comparison of (a) the MDCs at 0.12 eV above the DP, (b) the ARPES intensities around the DP, and (c) the EDCs 
at the P point between x = 1.0 and 0.9. The peak width of the MDC (Afc) obtained by the fitting with Lorentzian peaks is 
also indicated in (a). 
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FIG. S4: (a) Near-i?F ARPES intensity plots and (b) the corresponding EDCs for different material classes of topological 
insulators, Efo-^SbxTes-i/Sey [jiH and Pb(Bii_ x Sba;)Te4 Comparisons between the ordered (Bi2Te2Se, PbBi2Te4) and 

disordered (Bi1.5Sbo.5Te1.7Se1. 3, Pb(Bio.6Sbo.4)2Te4) phases are presented. The EDCs at the Y point are indicated by blue 
curves. None of them shows a Dirac gap. 



